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The orbital magnetic response of mesoscopic systems has been extensively studied in 
the last 25 years [I]. One of its manifestations is the occurrence of persistent currents 
and weak magnetism in small (sub-micron size) metallic rings threaded by a magnetic 
flux [2]. By 'small' one means that the circumference L of the ring is much smaller than 
the phase coherence length L^, and so quantum coherence is maintained throughout. 
This condition can be achieved at low enough temperatures. In a metallic ring disorder 
is very weak (as expressed in terms of the Ioffe-Regel condition kpl ^> 1, with kp 
and I being the Fermi momentum and the mean free path), and the currents persist for 
systems of a few microns in size [3], while for kpl < 1 they decay exponentially with the 
system size. The magnetic response of metallic rings proved to be an important tool to 
study fundamental aspects of quantum mechanics, such as quantum coherence [I] and 
the Aharonov-Bohm (AB) effect [HIE]. In systems exhibiting the AB effect, the magnetic 
flux appears as a pure Abelian U(l) gauge, in the sense that it only affects the phase of 
the wavefunction. Every observable quantity is therefore periodic in the magnetic flux 
(this is the Byers-Yang theorem [7]). Analogously, in the presence of a strong electric 
field (either internal or external), spin-orbit interactions give rise to the Aharonov- 
Casher (AC) effect [8] [9], which can be described as a pure non- Abelian SU(2) gauge. 




Figure 1. The sample considered in this work consists of two clean wires modeled as 
ideal rings with lengths Li and L 2 , areas A\ and A 2 , threaded by fluxes $! and $27 
touching at a contact point P. 

In this work, motivated by these issues, we examine the relevance of the topology 
of the sample to its magnetic response. We investigate thoroughly a closed system 
composed of two (ideal) metallic rings sharing a single contact point P (see Figured]). 
From a topological viewpoint, this sample has genus two, whereas a single ring has 
genus one. This difference turns out to enrich the physics in a rather non-trivial way, as 
already noticed in [10J. The present work includes a systematic study of the energy 
spectrum, persistent currents and magnetization of the sample shown in Figure (TJ 
including their dependence on the ring lengths, on the magnetic fluxes, and on the 
number N of electrons at zero temperature. The difference between the single and 
double- ring geometry is underlined throughout. The main novel features of the latter 
case are as follows. Wavef unctions generically depend on the fluxes. A rich pattern of 
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twofold and threefold degeneracies is observed. As a consequence, the question whether 
a given level has a paramagnetic or diamagnetic response is more subtle than in the 
single-ring geometry. The commensurability of both ring lengths is a key issue. In the 
case of commensurate ring lengths, a finite fraction of the levels are 'idle', in the sense 
that their energies do not depend on the magnetic fluxes, so that these levels do not 
contribute to the persistent currents. The case of incommensurate ring lengths requires 
a special treatment, inspired from the theory of modulated incommensurate structures. 
As for the AC effect, the two-ring geometry allows one to study a novel feature which is 
absent in the single-ring geometry. Suppose that the AC effect is realized by threading a 
ring with a long straight wire with constant longitudinal charge density. In a single-ring 
geometry this construction implies that s z , the spin component along the ring axis, is 
conserved, and so the problem decomposes into two independent ones, one for spin up 
and one for spin down. In the two-ring geometry, we can realize the AC effect as a pure 
gauge in an s z non-conserving system. Remarkably enough, this is possible only if the 
magnetic fluxes through both rings are non-trivial, i.e., neither integer nor half-integer 
multiples of the flux quantum. This kind of interplay between the AB and AC effects 
had not been noticed so far, to the best of our knowledge. Finally, the introduction of a 
non-Abelian SU(2) flux can also affect the sign of the sample's magnetization, turning 
a diamagnetic to a paramagnetic response or vice versa. 

In the following we elaborate and substantiate the issues presented above. Star- 
ting with the AB effect, the following topics are successively covered (section num- 
bers in parentheses): the Hamiltonian and its characteristic equation (j2J), basic 
observables, including persistent currents and magnetization various special cases of 
interest (jlj), twofold and threefold degeneracies (jHJ), and the spectrum and observables for 
commensurate (jSJ) and incommensurate ([7]) ring lengths. Then in Section [8] we examine 
the role of spin-orbit interactions and construct a Hamiltonian in terms of SU(2) fluxes 
leading to the AC effect. The energy levels and the magnetization are calculated in the 
special case of two equal rings, the emphasis being put on a novel AB-AC interference 
effect. A summary of our findings and a discussion are presented in Section [9l while two 
appendices are devoted to a reminder of the well-known case of a single ring (A), and 
to an extension of the analysis to three coupled rings (B). 

2. The Hamiltonian and its characteristic equation 

We consider a clean metallic sample in the form of two unequal rings touching at a 
contact point P, as shown in Figure HJ The rings are planar, but may otherwise assume 
arbitrary shapes. The rings have lengths L\ and L 2 and areas A\ and A 2 . In the presence 
of a uniform transverse magnetic field B, they are therefore threaded by magnetic fluxes 
$! = BAi and $ 2 — BA 2 . In the case of circular rings, to be used in numerical 
illustrations of our results, we have A\ = L\/(Aii) and A 2 = L|/(47r). In order to 
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compare both ring lengths, we introduce a variable < < 1 so that 
L\ L2 



(2.1) 



h\ + L 2 L\ + L 2 

It will turn out that the system has different characteristics for commensurate lengths 
(rational 0) and for incommensurate lengths (irrational 0). 

To write down the Hamiltonian, it is useful to employ reduced units (h = c = e = 
2m = 1), so that the flux quantum reads <3> — 27r. The spectrum and the observables 
are therefore 27r-periodic in $1 and <3? 2 . We parametrize a point of the left ring by 
its curvilinear abscissa < Si < Li, starting from the contact point P and oriented 
clockwise, and similarly a point of the right ring by < S2 < L2. Furthermore, we 
neglect spin degrees of freedom (except in Section [8]). The one-body Hamiltonian of the 
system reads 

n = { Pl - ai) 2 + (pa - a 2 ) 2 , (2.2) 

with 

Pi = -1 4-, P2 = -1 -j—, (2.3) 
as 1 as 2 

whereas the tangential vector potentials a\ and a 2 can be taken equal to 

ai = — , a 2 = — . (2.4) 

A state is described by a pair of wavefunctions {ip^(si), ip^ (S2)} ■ The first term 
of the Hamiltonian 7i acts on the left component ip^(si), whereas the second term 
acts on the right component ^ 2 )(s 2 )- The behavior of the wavefunctions at the contact 
point P is in general described by a unitary junction S-matrix. In the present work 
we make the simplest choice, which corresponds to just requiring the continuity of the 
wavef unction: 

^«(0) = V (1) (£i) = ^ (2) (0) = ^ 2) (L 2 ) = V>(P), (2.5) 
and the conservation of the current: 

(pi-ai)^ (1) (0)-(pi-ai)^ (1) M 
+ (Pa - a 2 )^ 2 \0) - (p 2 - a 2 )^ 2) (L 2 ) = 0. (2.6) 

Setting E = q 2 for the energy eigenvalue, we look for an eigenstate of 7i in the form 

^^(si) = e iaiSl (A 1 e igSl + ^e- 1 " 51 ), 

^( 2 )( S2 ) = e 102S2 (A 2 e 1<ZS2 + B 2 e- iqS2 ). ' ' ^ 

The condition (I2.5P allows one to express these four amplitudes in terms of ip(P) as 



1 2isingLx ^ ' 1 2i singly ^ ' 

-1*2 _ „-ig^2 p l 9^2 _ 0-1*2 

zi sin gL 2 2i sin gL 2 



(2- 



The condition (12. 6p then yields the characteristic equation 

D(q) = 0, (2.9) 
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where the characteristic function reads 

D(q) = sin q{L\ + L 2 ) — cos $2 sin qL\ — cos $1 sin qL 2 . (2-10) 

The eigenstates of the Hamiltonian correspond to the solutions q n > of ( 12. 9ft . 
labeled by n = 1,2,... and ordered as < q\ < q 2 < . . . The corresponding energy 
eigenvalues are E n — q\. In the special case where each ring is threaded by a quarter 
flux unit ($1 = $2 = 7r /2), we have q n = nn/{Li + L 2 ) (see (I4.10p ). In the general case, 
i.e., for arbitrary values of the fluxes, we set 

nir + q„ , . 

^ = T-rr (" = i,2..-0- (2.11) 

-Li + Li 2 

We will refer to g n as the modulation of the spectrum of eigenmomenta with respect to 
the linear behavior (14.101) . This quantity will be shown in Section [7] to obey the bound 

\g n \<7i. (2.12) 



The normalization of the wavefunction of the nth level reads 







Mi>n) = I ^\s 1 )-ds 1 + I ^(s 2 )-ds 2 = \MP)\ Z A( gn ),(2.13) 



2 



L-> 







(2), 



2 



i2 



where 



hence 



1 - cos$i cosg n Li 1 - cos $ 2 cosq n L 2 . , 

A(g„)=Li = + £2 = , (2.14) 

sin q n L x sin g n L 2 



^n(P) = A(g n )~ 1/2 , (2.15) 

up to an irrelevant phase factor. 

It is worth underlining a key difference between the present situation and that of 
a single ring, recalled in Appendix A. The gauge transformation employed in ( 12. 7ft . 
while it locally eliminates the vector potential from the Schrodinger equation, does not 
lead to a Bloch-type boundary condition, at variance with (IA.2I) . As a consequence, 
and in contrast with the single-ring wavefunction ip(s), the functions ^i(s) and ip 2 {s) 
given in (12.71) are not periodic in their respective arguments S\ and s 2 . The two-ring 
topology is indeed in marked difference with the single-ring one. Waves propagating in 
each ring are scattered at each passage at the contact point P. This multiple-scattering 
phenomenon destroys the periodicity of the plane waves characteristic of the single-ring 
problem. As a consequence, stationary states, as given by solutions of ( 12.91) . bear in 
general a non-trivial dependence on both ring lengths Li, L 2 and on both fluxes $1, $2- 
The periodicity of physical observables in the fluxes is however guaranteed by the Byers- 
Yang theorem, whose validity is independent of whether there is a Bloch analogue or not. 

3. Persistent currents and magnetization 

The contributions of the nth level to the persistent currents in each ring and to the 
magnetization read 

'* = -ip = ^ = - S -^ = A lh , n + A 2ha . (3.1) 
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The persistent currents can be evaluated as follows. Considering J l n for definiteness, 
we have 

T 1 9E » 2 / / I I / \ 

h,n = - ~r "5 — = T~ \Vn\Pi ~ ai\ip n ) 



^ jf'tf^i) (^d^" 01 ) ^^ds,. 



(3.2) 



Using the expression (12.71) of the normalized wavefunction, together with (12.81) . (12.141) 
and (I2.15p . we obtain after some algebra 

T -^1 n n T -^ n n A/T — O ( AlQl,n . A 2 Q2,n\ , qq n 

I l,n — ~T~ Vl,n; J 2,n — "T - V2,n; 1V1 n ~ Z Qn I j 1 j I i \<3.C>) 

L\ L 2 \ Li\ L 2 / 

where the dimensionless current amplitudes read 

_ Li sin$i _ L 2 sin $ 2 

Vi,n : j — 77 — r, V2,n : p— r~, — r- yoA) 

sm g n Li A(g n ) sm g n L 2 A(g n ) 

An alternative approach consists in evaluating the currents from the spectrum. We have 

dq n ( dD/d$>i\ sin$! sing n L 2 

h - = - 2 "- = 2q " UwL = 2tn wm, ■ (3 - 5) 



lq=q„ 



The current amplitudes can therefore be expressed as 

n t &Qn dg n dq n dg n 

Q lt n = -Li — — = -0— — , Q 2 .n = -^2 777- = "(I - 0) (3-6) 

The following identity ensures that the results (13.31) . (13 .4p and (13. 5p . (13. 6p are identical: 
(ir) = -sing n L! sing n L 2 A(g n ). (3.7) 

\ y / q=q n 

For a zero-temperature system with N electrons, the N lowest energy states are 
occupied. The total energy and magnetization are therefore given by 

N op N 

E="£ E n, M = -7^ = H M n- (3-8) 
n=l dB n=l 

Finally, explicit bounds on various quantities of interest can be derived by applying 
the inequality | (ij; n \0\ij; n ) | 2 < (ij} n \0 2 \ij} n ) to the operators 0i — p\ — a\ and 2 = £> 2 — a 2 . 
The expression (13. 2p implies that the persistent currents obey the bound 

(L 1 / 1 , n ) 2 +(L 2 7 2 , n ) 2 <4^ ) (3.9) 

i.e., 

Qln + Qln < 1- (3-10) 

Furthermore, using (13.51) and (13. 7p . we can respectively recast (13. 9p as 

^>(^^)+(^P)' (3.11) 
\ sm q n Li ) \ sm q n L 2 J 

and 

/ dD \ 2 

I — ] > (L x sin $! sin q n L 2 ) 2 + (L 2 sin $ 2 sin q n L^ 2 . (3.12) 



<9g 



q=q n 
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4. Special cases of interest 

4-1. No magnetic fluxes 

We consider first the problem in the absence of fluxes: $1 = $2 = (mod. 27r)JJ| The 
characteristic function f l2.10p factors as 

r>r \ a ■ <?(£i + £2) . qL x . qL 2 

V(q) = —4 sin sin — — sm — — . (4.1) 

The spectrum therefore consists of the following three sectors, in correspondence with 
the factors of the above expression. 

• Bilateral states. These states correspond to sing(Lx + L 2 )/2 = 0, hence 

q = T^rr (6 = 0,2,4,...), (4.2) 

Li\ -+- Li 2 

with '6' for bilateral. The corresponding wavefunctions are standing waves living 

on the whole system: 

id) ( \ /m s cosgsi + cosg(Li - si) 

^W(si) = V(P) r- 7 , 

1 + cosgLi ^ 43 n 

1(2)1 \ , m , cos qs 2 + cos q(L 2 - s 2 ) 

T 2) {s 2 ) = V(P) 7- j • 

1 + cosqLi 

• Left states. These states correspond to singLi/2 = 0, hence 

q= 1 ^- (/ = 2,4,6,...), (4.4) 

with T for left. The corresponding wavefunctions are standing waves living on the 
left ring, whose amplitude vanishes at the contact point § 

ip {1 \s 1 ) ~ singsi, tp {2) (s 2 ) = 0. (4.5) 

• Right states. These states correspond to singL 2 /2 = 0, hence 

rix 

q=— (r = 2,4,6,...), (4.6) 

with V for right. The corresponding wavefunctions are standing waves living on 
the right ring, whose amplitude vanishes at the contact point: 

^ {1) (si) = 0, ^ 2 \s 2 ) ~ sings 2 . (4.7) 



The ground-state energy, obtained by setting b = in ( 14.21) . vanishesljjj The 
modulation introduced in (12. 1 1 P reads g± = —it, which saturates the bound (I2.12p . 
The corresponding ground-state wavefunction is uniform: ip (si) = ip ($2) = 

For small fluxes, we have 

| We recall that the notation (mod. 2tt) means up to a multiple of 2tt. 

§ The symbol ~ is used whenever the wavefunction normalization is not given explicitly. 

II Notice that I = and r = are not allowed in (HH) and (l4~6l). 
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so that the ground state is always diamagnetic. As far as excited states are concerned, 
the results (14. 2p . (I4.4p and H4.6[) respectively become 

$2 _ $2 

q(Li + L 2 ) ^ for + — ^ — - cot 6tt0, 

qL x ^ Zvr + cot — , (4.9) 

2 

qLo ~ r7r H cot . 

y 2 1-0 

These results show that there is no general rule to predict whether a given (left or 
right) state is paramagnetic or diamagnetic. Bilateral states are always hybrid (one 
paramagnetic ring and one diamagnetic one). 

The only degeneracies between the three interlaced spectra f)4.2j) . (14.41) and (14.61) 
are the threefold ones taking place in the commensurate case at qa = 9 = ttu, where 
v = 2yU = 2, 4, 6, . . . is an even integer. With the notations of Section [61 these momentum 
values correspond to b = fim, I = jumi and r = fim 2 . The quadratic terms in (14.91) 
diverge. This is in agreement with the fact that near a degeneracy the energy eigenvalues 
vary linearly with the fluxes, rather than quadratically (see Section [5] for more details). 

4-2. Integer or half-integer fluxes 

Whenever each flux is either integer or half- integer (i.e., an integer or a half-integer 
multiple of the flux quantum <3> = 27r), the spectrum of the system still consists of the 
above three types of states (bilateral, left and right). The corresponding momenta are 
still given by (I4.2p . (14.41) and (14.61) . Table[T]lists the characteristics of the spectrum in the 
four different cases, corresponding to $i = or it (mod. 2ir) and $2 = or it (mod. 2tt). 



$1 


$ 2 


D(q) 


b 


/ 


r 








. q(Li + L 2 ) qL x . qL 2 

—4 sin sin sin 

2 2 2 


0,2,4,... 


2,4,6,... 


2,4,6,... 





7T 


q{Li + La) . qL x qL 2 

4 cos sin cos 

2 2 2 


1,3,5,... 


2,4,6,... 


1,3,5,... 


7T 





q(Lx + L 2 ) qLx qL 2 

4 cos cos sin 

2 2 2 


1,3,5,... 


1,3,5,... 


2,4,6,... 


7T 


7T 


. q(Lx + L 2 ) qLx qL 2 

4 sin cos cos 

2 2 2 


2,4,6,... 


1,3,5,... 


1,3,5,... 



Table 1. Characteristics of the system for the four cases with integer or half- integer 
fluxes: factorized form of the characteristic function D(q); integers b, I and r entering 
the spectra (|4.2j) , (|4.4[) and ()4.6j) , respectively corresponding to bilateral, left and right 
states. 



Left and right states are met in a more general setting. More precisely, left 
states with even (resp. odd) / exist as soon as $1 = (resp. $1 = it (mod. 2tt)), 
irrespective of $2; whereas right states with even (resp. odd) r exist as soon as $2 — 
(resp. $ 2 = 7T (mod. 2ir)), irrespective of $1. 
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4-3. Quarter-integer fluxes 

The situation where each ring is threaded by a quarter flux unit ($1 = $ 2 = 7r/2) is 
also a special case of interest. In this case, the characteristic function boils down to 
D(q) = sinq(Li + L 2 ). The momenta therefore have the linear dependence 

Q" = T— T (n = l,2,...). (4-10) 

In spite of the simplicity of the spectrum, the current amplitudes have the following 
non-trivial expressions: 

Qi,n = -0sinwr0, Q 2 ,n = (-l) n (l - 0) sinn7T0. (4.11) 

These quantities obey \Qi >n \ < 0, |Q2,n| < 1 — 0, and \Qi >n \ + lQ2.nl < 1; the latter 
bound being more stringent than the general one ( I3.10p . The currents in both rings 
therefore have the same sign (resp. opposite signs) whenever the level number n is odd 
(resp. even). 

Using (13.61) . the result (I4.1ip can be recast as follows. If the fluxes are close to a 
quarter flux quantum, setting $1 = n/2 + <5$i, $ 2 = tt/2 + 5$ 2 , the modulation g n 
introduced in ( 12. lip reads, to first order in 5$i and £$ 2 , 

f 5&i — 5&2 for n even, , 1 _ nS 

o n Rismn7r0x (4.12) 
\ <5$i + <5$ 2 for n odd. v y 

5. Degeneracies 

A degeneracy manifests itself as a multiple (i.e., at least double) root of the characteristic 
equation ( 12. 9p . i.e., as a simultaneous root of D(q) = and dD/dq = 0. The 
bound (I3.12p shows that the system has no accidental degeneracy. A degeneracy may 
indeed only take place when both products sin $1 sin gL 2 and sin $ 2 sin qL\ vanish 
simultaneously, i.e., when at least one factor of each product vanishes. It can be checked 
that there are only two types of degeneracies, to be successively investigated hereafter. 

5.1. Twofold degeneracies 

Twofold degeneracies occur when either sinqLi = sin$x = (but singL 2 7^ 0) or 
singL 2 = sin$ 2 = (but singLi 7^ 0). These two instances will be respectively referred 
to as left and right degeneracies, for a reason that will become clear in a while. 

Consider the first instance for definiteness. This is a left degeneracy, as the state 
to become twofold degenerate is a left state, in the sense of Section 14.11 We have 
cos (/Li = cos$i = (—1)', with the notation (14.41) . so that 
3D 

— = (-l)'Li(cosgL 2 - cos$ 2 ). (5.1) 
oq 

The conditions for a twofold degeneracy are therefore 

sin qL\ = sin $1 = 0, cos qL 2 = cos $ 2 . (5.2) 
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At the degeneracy the momentum is q = Ik/Li, whereas the fluxes read $x = (2j + 1)tt 
and $2 = {2k+elL 2 / Li)n, where e = ±1, whereas j and k are integers. For circular rings 
(more generally, for rings with similar shapes), we have $2/^1 = A2/A1 = [L 2 jL\) 2 = 
((1 — 0)/0) 2 . This condition relates to e, j and k as follows: 

e = +l: 2 (j-A; + /)0 2 -(4j + 3O0 + 2j + / = O, 

£ = -1 : 2(j - k)0 2 - (4j + /)0 + 2j + Z = 0. 1 ' ' 

The system may therefore have twofold degeneracies both in the commensurate case (0 
rational) and in the incommensurate case provided is a quadratic number, obeying an 
equation of the form 

I0 2 + J0 + K = 0, (5.4) 

where /, J and K are integers. 

In order to explore the vicinity of the degeneracy, we set 

g= — +?7, $i = /vr + 5$i, $2 = -^ h<5$ 2 - (5.5) 

By expanding the characteristic equation (12. 9p to second order in 77, <5$x and 5§2, we 
obtain the reduced equation 

+ 2L 2 )r/ 2 - 2L 1 5$ 2 r/ - <5$ 2 = 0. (5.6) 

The momenta of the two nearly degenerate states are therefore given by 

with the definition 

W 2 = (Lx (Li + 2L 2 )<5$ 2 + L 2 5$ 2 ) 1/2 . (5.8) 

The momenta shifts rj± vanish linearly with the distance to the degeneracy in the 
(6&i,6<&2) plane. The corresponding amplitudes of the persistent currents read 

^ L 1( 5$i L 2 / L 1( 5<l>2\ n . 

These amplitudes therefore remain of order unity, and they depend on the direction in 
the (<5$i,<5$2) plane along which the degeneracy point is approached. As this direction 
is varied, the amplitudes describe an ellipse: 

(1 - 0) 2 Ql + 0(2 - 0)Q\ - 20(1 - 0)Q 2 = 0. (5.10) 
5.2. Threefold degeneracies 

We have noticed at the end of Section H~T1 that the spectrum has threefold degeneracies in 
the absence of fluxes in the commensurate case at qa = 9 = iru, where v = 2/x = 2, 4, . . . 
is an even integer, with the notations of Section [61 
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More generally, there are threefold degeneracies at all the energies of the idle states, 
given by (I6.4p . i.e., qa = = ttu, with v = 1,2,... Consider indeed one of these energies. 
We have 

3D 

— = ae 1 e 2 (m 1 (l - e 2 cos$ 2 ) + m 2 (l - e x cos$i)), (5.11) 
with the notations (16. 6p . The conditions for a threefold degeneracy are therefore 

COS$i=£i, COS $2 = ^2- (5.12) 

Right at the degeneracy, the fluxes therefore read $i = numi and $2 = Ttvm 2 (mod. 2tt). 
The situation in the absence of fluxes is recovered when v is even. If v is odd, at least 
one of the fluxes is a half integer, because at least one of the integers mi, m 2 is odd. 
In order to explore the vicinity of the degeneracy, we set 

6 = vii + m $1 = Tivmi + 6<&i, $2 = ixvm 2 + <5$ 2 - (5.13) 

By expanding the characteristic equation ( 12. 9ft to third order in 77, and to second order 
in £$1 and 8<&2, we find three solutions: the idle state at rj = 0, and two symmetrically 
shifted active states at r)± = ±W3, with 

fwg+miwiy* (514) 

\m 1 m 2 {mi + m 2 ) J 

The momentum shift W3 again vanishes linearly with the distance to the degeneracy in 
the £$2) plane. The corresponding amplitudes of the persistent currents read 

Qi,± = =F , , \ w , Q2,± = T ( £ 2 (5.15) 

Here again, these amplitudes remain of order unity and depend on the direction in the 
(8^1,8^2) plane along which the degeneracy point is approached. As this direction is 
varied, the amplitudes describe an ellipse: 

(1 - 0)Q\ + 0Q 2 2 - 0(1 - 0) = 0. (5.16) 



6. Commensurate ring lengths 

We now turn to the situation where the ring lengths L\ and L 2 are commensurate. 
One of the peculiar features of this situation is the existence of threefold degeneracies, 
already investigated in Section 15.21 

In the commensurate situation, both ring lengths are multiples of the same 
fundamental length a: 

Li = mia, L 2 = m 2 a, (6.1) 

where the integers mi > 1 and m 2 > 1 are relatively prime. Setting m = m\ + m 2 > 2, 
the variable takes the rational value = rai/m. For circular commensurate rings, the 
fluxes read $1 = ml Ba 2 / (47r), $2 = fn 2 Ba 2 / (47r), so that the magnetization is periodic 
in the magnetic field B, with period B = 8ir 2 /a 2 . 
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6.1. Spectrum 

Introducing the reduced momentum 

9 = qa, (6.2) 
the characteristic function becomes 

D(9) = smm9 — cos $2 sinmi# — cos$i sinm2#. (6-3) 
The spectrum of the system is 27r-periodic in 9. It consists of two types of states. 

• Idle states. They correspond to the trivial solutions of the characteristic equation: 

Q = txv {v = 1,2,...). (6.4) 

These states carry no current, as their energy is independent of the fluxes. The 
corresponding wavefunctions do however depend on the fluxes. We have indeed 
(si) ~ (e- 1 * 2 / 2 - e 2 e 1$2 / 2 )e iaiSl sin qs u 
¥ 2) (s 2 ) ~ fci e '* l/2 - e-^/ 2 )e m ^ sin qs 2 , ( ' 



with 

e x = e iqLl = {-l) vm \ e 2 = e iqL2 = (-iy m \ (6.6) 

Idle states become threefold degenerate for integer or half-integer values of the 
fluxes such that the condition (15.121) is fulfilled. 

• Active states. These states, corresponding to the other solutions of the character- 
istic equation D{9) = 0, carry non-zero persistent currents in general. 
Setting c = cos^, we have the identity sin m6 = sind ?7 m _i(c), where U n (c) is the 
nth Tchebyshev polynomial of the second kind, whose degree is n. Focusing onto 
active states, the characteristic equation therefore reads 

f/ m _i(c) - cos $ 2 t/ mi _i(c) - cos $1 t/ m2 _i(c) = 0. (6.7) 

This equation has m — 1 solutions C(k) {k = 1, . . . , m — 1). We write these solutions 
as C(fe) = cos#(fc), with < 9qa < ir, ordered as 

< 9 {1) <...< 9 {m . x) < n. (6.8) 

Threefold degeneracies correspond to the limiting situations 9^ = or = tt, 

whereas twofold ones can take place anywhere along the sequence of 0(fc)'s. 

To sum up, in each period of the spectrum of length 2ir in the variable 9 = qa, there 
are 2(m — 1) active levels and 2 idle ones, i.e., a total of 2m states. The modulation and 
the current amplitudes associated with the states at 9 = 9^) (k — 1, . . . , m — 1) read 

g k = m9( k) - kir, Q hk = -mi Q2,k = ~ m 2 (6.9) 

whereas those associated with the states at 9 = 2n — 9n\, 

g2m-k = kn - m9(k), Qi^m-k = mi Q 2 ,2m-k = m 2 ^p", (6-10) 
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are the opposites of the first expressions. The bound (12.121) implies 

(k - 1)tt < m9 (k) < (k + 1)tt (k = 1, . . . , m - 1). (6.11) 

Modulation and current amplitudes then repeat themselves periodically, with period 2m. 
The main characteristics of the states in the first period are listed in Table [2j 



n 


1 




m — 1 


m 


m + 1 




2m -1 


2m 


6 


0(1) 




9(m-l) 


IT 


27T — 0( m _i) 




2ir - (1) 


2ix 


type 


A 




A 


I 


A 




A 


I 


9n 


9i 




9m-l 





—gm-i 









Ql,n 


Qi,i 




Ql,m—1 





— Ql,m-1 




-<5i,i 





Q2,n 


Q%i 




Q2,m-1 





—Q2,m-l 




-Q 2 ,i 






Table 2. The first 2m states of the system in the commensurate case, corresponding to 
the first period of the spectrum in the variable 9 = qa. For each column, corresponding 
to a state, the table gives the level number n, the angle 9, the type of state (T for idle 
or 'A' for active), the modulation and the current amplitudes. 



6.2. Rings with equal lengths 

The case where both rings have equal lengths is the simplest of all the commensurate 
situations. We have m\ = m2 = 1, m = 2, = 1/2 and a = L = L\ = L2. The 
characteristic equation (16.71) has a single solution: 

cos $! + cos $ 2 

c ( i) = cos = . (6.12) 

Idle states (n = 2p) and active states (n = 2p — 1) alternate along the spectrum: 

W (2 P -l)7T + (-l)P- 1 ff 1 

<? 2p = gap-i = 2£ ' ( 6 - 13 ) 

The modulation and the current amplitudes of the lowest active state (n — 2p — 1 = 1, 
i.e., p = 1) read 

^ sin $1 ^ sin $ 2 

In the case of two equal circular rings, the above predictions can be made fully 
explicit. We have $1 = $2 = $ = BL 2 /(4tt). Symmetry considerations allow us to 
restrict ourselves to < $ < n, so that 0(i) = $, ^1 = 2$ — 7r and Q11 = Q 2> i = —1/2. 
Hence, using (I3.3p . the contribution of any odd state (n = 2p— 1) to the magnetization 
reads 

M 2p _ 1 = (-ir^ + ^. (6.15) 

Inserting this formula into (13.81) . we obtain the expression of the magnetization of a 
system with electrons at zero temperature. This expression, given in Tabled depends 
on A^ (mod. 4). The corresponding result for a single circular ring, given by (IA.6[) . (I A. T[) . 
is also listed in the Table for comparison (the latter depends on the sign of N, i.e., 
on N (mod. 2)). 
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N (mod. 4) 


M (two equal rings) 


M (one single ring) 





N(tt -$)/(4tt) 


N(n - $)/(2tt) 


1 


-(iV + 1)$/(4tt) 


-JV$/(27r) 


2 


-JV$/(47r) 


A^(tt - $)/(2?r) 


3 


(A^ + l)(7T-$)/(47T) 


-JV$/(27r) 



Table 3. Expression of the magnetization M of a system with TV electrons at zero 
temperature, as a function of the flux <I> in the range < $ < n, for two equal circular 
rings and for a single circular ring. 



6.3. Magnetization 

The simple linear dependence of the magnetization on the magnetic field in the range 
< $ < 7r, observed in Table [31 is a peculiarity of the cases considered there, namely 
one single circular ring or two equal ones. 

In all the other cases of commensurate circular rings, the curve M(B) is non- 
trivial. This is illustrated in Figure EJ showing plots of the magnetization over one 
period, against $i/(2tt) = B/B = BLl/(8ir 2 ), for L 2 = L x (left) and L 2 = 2L X (right). 
The latter example (mi = 1, ni2 = 2, m = 3) is generic of the commensurate case, 
except that the characteristic equation (I6.7P is of degree two in c and can therefore be 
solved analytically. We thus obtain 

= - (cos $i + (4 + 4 cos $ 2 + cos 2 $i) 1/2 ) , 
= - (cos $1 - (4 + 4 cos $ 2 + cos 2 $i) 1/2 ) . 



C(l) 



C(2) 



COS (7(1) 
COS 01 



(6.16) 



The currents and the magnetization then follow from (16.91) . (I6.10p and (13.31) . (13.81) . 





0.4 0.6 
* /(2tt) 



Figure 2. Plots of the magnetization M over one period, against $i/(27r) = 
BLf/(8TT 2 ), for commensurate circular rings with L 2 = Li (left) and L 2 = 2L X (right), 
for N — 3 (empty circles), N — 4 (full circles) and N = 5 (full triangles). Data for 
N = 3 and N = 4 coincide on the left plot, as n = 4 is an idle state. 
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7. Incommensurate ring lengths 

We now turn to the situation where the ring lengths L\ and L 2 are incommensurate, 
i.e., the variable is irrational. 

The key quantity is again the modulation g n of the spectrum. In the commensurate 
case, i.e., for a rational = mi/m, g n has been shown to be periodic in n, with period 2m. 
In the present incommensurate case, g n is expected to never repeat itself exactly. At a 
quantitative level, the structure of g n is revealed by the situation where $1 = $2 = 7r/2, 
considered in Section 14.31 The result (I4.12p suggests the Ansatz 

n even : g n = G c (nn0), 
n odd : g n = G o (n7T0), 

where G c ('e' for even) and G Q ('o' for odd) are 27r-periodic functions of x = nir0. The 
modulation is therefore a quasiperiodic function of the level number n. We will refer 
to G c and G Q as the hull functions, following the term introduced by Aubry in the 
context of modulated incommensurate structures [H] . The current amplitudes are then 
also quasiperiodic functions of n: 



dG e (nir0) /-in dG e (nm 

neven: Q l n = -0 — , Q 2 ,„ = -(1 - 0) — — 

<9G o (mr0) dG (nvr£ 

nodd: Q 1>n = -0 — , Q 2 ,n = -(l-0) 



(7.2) 



Inserting (17. II) into the characteristic equation (12.91) leads to implicit equations for 
the hull functions: 

sin G c = cos $ 2 sin(x + 0G e ) — cos $1 sin(x — (1 — 0)G e ), . . 

sin G Q = — cos $ 2 sin(x + 0G O ) — cos $1 sin(x — (1 — 0)G O ). 

These equations imply the following properties: G e (x) and G Q (x) are 27r-periodic, odd 
and continuous functions of x\ G c (x) = for cos$i = cos $2; whereas G Q (x) = for 
cos$i = — cos $2; changing cos $2 into its opposite amounts to exchanging G c and G Q , 
whereas changing cos $1 and cos $2 into their opposites amounts to changing x into x+n. 



7.1. Integer or half -integer fluxes 

It is worth investigating first the case where there are no magnetic fluxes, and more 
generally the situation where the fluxes are integer or half-integer. The spectrum in the 
absence of magnetic fluxes has been studied in Section 14.11 

• Even values n = 2p of the level number correspond to bilateral states. The 
expression (14. 2 p shows that the modulation vanishes whenever n = b = 2p is an 
even integer. We have therefore G c = 0. 

• Odd values n = 2p — 1 of the level number correspond to left and right states. For 
left states, setting / = 2/, ( T4~4l) yields 2(p-l)n < (2p - l)n + g 2p -i = 2ln/0 < 2pn, 
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(7.4) 



Similarly, for right states, setting r = 2r, ( 14. 6ft yields 2(p—l)ii < (2p— l)7r+^2p-i 
2r7r/ (1 — 0) < 2p7r, hence 



P 



1 + Int 



5f 2p -i = 2vr Frac 



(7.5) 



1-0 V 1-0 

The expressions ( 17. 4p and (17.51) cover every integer p once. The inverse formulas 



read 



Fracp0 < : I — Intp0, 



92p-l 



Fracp0 >0: r = p — 1 — Intp0, g 2p - 



71 



7T — 



2ir Fracp0 

' 
27r(l — Fracp0) 

1-0 ' 



(7.6) 



The latter expressions for the modulation can be brought to the form (17.11) with 
G Q (x) = Q(x), the periodic, odd, piecewise linear continuous function defined for 
< x < 2tt as: 

{— x/0 for < X < 7T0, 

(x-7T)/(l-0) for 7T0 < X < tt(2 - 0), (7.7) 
(2vr - x) /0 for (2 - 0)vr < x < 2tt. 

The linearly increasing (resp. decreasing) parts of Q{x) describe right (resp. left) states. 
The cusps at x = tt0 and x = tt(2 — 0) eventually correspond to threefold degeneracies. 




-1.0 



0.2 0.4 0.6 0.8 1.0 
x/(27r) 

Figure 3. Plot of the hull functions G(x) (full line) and Q(x + ti) (dashed line), divided 
by 7T, against x/(2w) over one period, for = 1/t 2 . 



More generally, for the four situations corresponding either to integer or half-integer 
fluxes, considered in Section H~2| the non-zero hull functions are either equal to Q{x) or 
to Q{x + it), as shown in Table HI Figure [3] shows a plot of these two functions in the 
prototypical example of an incommensurate situation, namely L 2 /Li = r, where r = 
(l + V5)/2 w 1.618 034 is the Golden mean, so that = 1/r 2 = (3- V5)/2 w 0.381966. 

We recall that the integer part Int x and the fractional part Frac a; of a real number x arc defined by 
x = Into; + Fracx, with Int x integer and < Frac a; < 1, so that Fracx is periodic in x, with unit 
period. 
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$1 


$ 2 


G e (x) 


G (z) 











<?(*) 





TT 


£(x) 





TT 





+ 7r) 





TT 


TT 





Q(x + 7r) 



Table 4. Expression of the hull functions G c (x) and G (x), respectively characterizing 
the states with even and odd n, for the four cases with integer or half-integer fluxes. 
The function Q(x) is defined in (J7.7D . 



7.2. The general case 

Coming back to the general case, i.e., arbitrary values of the fluxes, we are now in a 
position to show another remarkable property: the hull functions are bounded by their 
limiting values for integer or half-integer fluxes: 

0<x<tt: Q{x) < G e (x), G Q (x) < Q{x + tt), 

tt<x<2tt: G(x + 7t) < G c (x), G Q (x) < Q{x). [ ' } 
In other words, the hull functions are inscribed in the two parallelograms shown in 
Figure [3j This property implies in particular |G e (x)| < tt and |G (x)| < tt, hence the 
bound (12321) . 

The inequalities ( jT.8|) can be proved as follows. It can be checked, using their 
expressions ( \7.2L that the current amplitudes Qi )Tl and Q2, n have well-defined signs: 

sign(/i in ) = sign(Qi >n ) = - sign(sin $i) sign(sinmr0), 

sign(7 2) n) = sign(Q 2 ,n) = - sign(sin $ 2 ) sign(sinri7r(l - 0)), (7.9) 

(— l) n— 1 sign(sin niz0) 

as long as both sin $1 and sin $2 are non-zero, so that degeneracies are avoided. These 
signs are constant over the domain < $1 < tt, < $2 < 7T. In particular, the 
observation made in Section 14.31 that the currents in both rings have the same sign 
(resp. opposite signs) whenever the level number n is odd (resp. even), holds all over 
this domain. The hull functions therefore have a monotonic dependence on $! and $2 
in the same domain, and they take their extremal values at the corners of the domain, 
i.e., for integer or half-integer fluxes. Figure H] shows plots of the hull function G Q for 
= 1/r 2 , both for $ a = $ 2 (left) and for $! = and variable $ 2 (right). 

If one flux is integer or half-integer, albeit the other is not, the hull functions G e (x) 
or G (x) exhibit linear parts, where they coincide either with G(x) or with Q[x + tt). 
These linear parts describe left or right states. They end at cusps which eventually 
correspond to twofold degeneracies. Consider for definiteness $1 = and < $ 2 < tt. 
The hull functions start linearly as: 

n even : G c (x) = —x/0 for < x < x c = 0$ 2 , ,-, in x 

n odd : G Q (x) — —x/0 for < x < x Q = 0{tt - $ 2 ). 

Let us close up this section with some numerical illustrations of our results in 
the case of circular rings. The observables (persistent currents and magnetization) are 
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x/(2tt) x/(2tt) 

Figure 4. Plots of the hull function G D (x), divided by n, against x/(2tt) over one 
period, for = 1/t 2 . Left panel: $i = <J> 2 = A;7r/10. Right panel: $i = 0, <& 2 = fc7r/10. 
The latter case exhibits the linear parts (|7.10|) . In both cases, k = 0, . . . , 10, bottom 
to top in the left part of the curves. 



given in terms of the hull functions by (17. 2L and (13. 3ft . (13.81) . Figure [5] shows plots of 
the current amplitudes Q\, n and Q2, n of individual levels against level number n, for a 
system of two circular rings with $i = £>L 2 /(47r) = ir/3, in a typical commensurate case 
(left): L 2 = 2Li, i.e., = 1/3, and in a typical incommensurate case (right): = 1/t 2 . 
The period 2m = 6 predicted in Section [6] is clearly observed in the commensurate case. 
Figure [6] shows plots of the magnetization against $i/(27r) = BLl/(8ir 2 ), for N = 3 
(left) and N = 10 (right), for two incommensurate cases corresponding to the nearby 
irrationals = 1/r 2 ~ 0.381966 and = 1/e ~ 0.367879. Twofold degeneracies are 
observed in the first case, in agreement with (15. 4p . as = 1/r 2 obeys 2 — 30 + 1 = 0. 




n n 



Figure 5. Plots of the current amplitudes Qi :U (full symbols) and Qi, n (empty 
symbols) against level number n, for 4>i = BL\j{Am) = tt/3. Left: commensurate 
case Z/2 = 2Li, i.e., = 1/3. Right: incommensurate case = 1/t 2 . 
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0.5 1.0 1.5 2.0 2.5 3.0 



-2 







0.5 1.0 1.5 2.0 2.5 3.0 



$/(2tt) 




Figure 6. Plots of the magnetization M against $i/(27r) = BL 2 /(87r 2 ), for N = 3 
(left) and N = 10 (right). Full lines: = 1/r 2 (jumps due to twofold degeneracies are 
shown by symbols). Dashed lines: = 1/e. 



8. The role of spin-orbit interaction 

So far spin did not play any role in our discussion. We now investigate the case that, 
besides the Abelian U(l) magnetic fluxes $i and $2; there are also non-Abelian SU(2) 
fluxes \&i and 

For the sake of consistency, let us recall some basic notions pertaining to the physics 
of SU(2) fluxes. Such fluxes arise as a result of spin-orbit interaction. Unlike U(l) 
fluxes, SU(2) fluxes are invariant under time reversal. While the U(l) flux leading to 
the Aharonov-Bohm effect is realized by threading a ring with a magnetic field, the 
SU(2) flux leading to the Aharonov-Casher effect is generated by piercing a ring with a 
line of charge. More precisely, if a system of electrons is confined to a plane and subject 
to an electric field generated by a straight perpendicular charged wire with constant 
charge A per unit length, we have the SU(2) analogue of the Aharonov-Bohm effect. 
The starting point of the analysis is the Pauli Hamiltonian. In the presence of a vector 
potential A and of an electric field E, within the approximate U(l)(g)SU(2) symmetry of 
the non-relativistic Schrodinger equation [12], this Hamiltonian reads, in dimensionful 
form, 



where 7 = e/(4mc 2 ). The third term on the right-hand side of (18 .ip is responsible 
for the spin-orbit interaction. In the case of a circular ring of radius R pierced by a 
charged wire through its center, we have E = 2Xr/R 2 , whereas the curvilinear abscissa 
and the circumference read s = R9 and L = 2ttR. In this geometry, the U(l) and 
SU(2) potentials appearing in ( 18.1 ft can be eliminated locally by the respective gauge 




(8.1) 
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transformations 

le r s \ ( 



, u(1) =exp(-£/Vdr)=exp 
9su(2) = exp ^-17 J (E x a) ■ dr^j = exp 



n\ ( 8 - 2 ) 



2ix J 1 

where the integrations are carried out along the ring. In the above expression for (?st/(2); 
the integral need not be path-ordered. We have indeed (E x a) ■ dr = — (E x dr) ■ a = 
-(2X/R 2 )(r x dr) - a = -2An ■ trdO, where n is the (properly oriented) unit vector 
perpendicular to the plane of the ring. As a result, the dimensionless SU(2) flux reads 
^ = — 47T7A. The above simplifying property is specific of the case where the electric field 
lies in the plane in which the electrons are confined. It has also been used in a spherical 
geometry, in deriving a tight-binding version of the L- S spin-orbit interaction [13]. The 
SU(2) flux appears as a 'pure gauge', in the sense that it only affects the phase of the 
wavefunction. However, the U(l) and SU(2) potentials cannot be eliminated globally. 
These potentials rather give rise to non-integrable phase factors, e 1 * for U(l) and e 1 * TV<T 
for SU(2). Observables are therefore 27r-periodic both in $ and in ty. The upshot of the 
above discussion is that, in the ring geometry, the Pauli Hamiltonian can be replaced 
by a simpler one, namely 

H = (p-a-bn- erf , (8.3) 

with the notations p = — id/ds, a = $/L and b = ty/L. For a single ring in the (x,y) 
plane, we have n = z, so that s z is conserved. The problem then reduces to that of two 
independent Aharonov-Bohm systems of polarized electrons, one with flux $ + $ and 
the second with flux $ — 

A somewhat less simple manifestation of the AC effect in terms of SU(2) fluxes 
occurs when s z is not conserved. This can be realized (for example) by considering 
two rings in different planes respectively subject to perpendicular magnetic fields and 
pierced by perpendicular lines of charges with uniform densities Ai and A 2 . Such a 
Gedankenexperiment is schematically displayed in Figure Thus, while in the U(l) 
two-ring problem the geometric configuration of the two rings (in particular their relative 
orientation) did not play an important role, the sample's geometry becomes of central 
importance when SU(2) fluxes are considered. This situation allows one to study the 
influence of SU(2) fluxes on the U(l) magnetization [14]. The main messages of the 
analysis given below are as follows, (i) Even when fix is antiparallel to ri-2, the non- 
conservation of s z requires both U(l) fluxes to be non-trivial, i.e., neither integer nor half- 
integer. This is a novel situation where the energy levels are sensitive to a combination 
of AB and AC effects, (ii) Even though ^ enters as a pure gauge, spin-orbit interaction 
might change the sign of the magnetization, so that the system switches between a 
paramagnetic and a diamagnetic response, (iii) In complete analogy with the U(l) 
case, where the orbital magnetization is related to the derivative of the ground-state 
energy with respect to $, it is natural to define an ! SU(2) magnetization' related to the 
derivative of the ground-state energy with respect to This magnetization is another 
equilibrium property. In each ring and for a given level n, it is proportional to the 
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Figure 7. Gcdankcncxpcrimental realization of the two rings with U(l) and SU(2) 
fluxes such that spin is not conserved. Solid lines represent magnetic fluxes $1 and 
3>2 as well as charged wires with respective longitudinal charge densities Ai and A2. 
Full arrows mark the direction of the electric fields, whereas empty arrows indicate the 
directions of the spins (more precisely, ni,2)- This is the special case corresponding to 
the Hamiltonian (|8.4p where the two axes ni^ lie in the same (x, z) plane. 



expectation value of the commutator {v, h ■ er}, where v is the velocity operator. Thus, 
it can in principle be measured, despite the fact that the SU(2) (spin) current is not 
conserved while the U(l) (charge) current is. 

The one-electron Hamiltonian of our two-ring system now reads 

H = (pi - ai - b\ -hi ■ er) 2 + (p 2 - a 2 - b 2 h 2 ■ cr) 2 , (8.4) 

where p\ and p 2 are the differential operators of (I2.3|) . hi and h 2 are two arbitrary unit 
vectors, cr is the vector of Pauli matrices, whereas the U(l) vector potentials a\ and a 2 
and their SU(2) analogues b\ and h 2 are related to the corresponding fluxes as follows: 

a 1 = — , a 2 = — , 61 = — , b 2 = — , (8.5) 

L\ L 2 L\ L 2 

The problem mostly depends on the angle (3 between both directions hi and h 2 in spin 
space, such that hi ■ h 2 = cos j3. For convenience we choose axes such that hi = (0, 0, 1) 
is along the z-axis whereas h 2 = (sin /3, 0, cos /3) is in the (x, z) plane. 

A state is described by a pair of wavefunctions {ift^(si), iJj^(s 2 )}, each of 
them being a two-component spinor. Separating spin and orbital degrees of freedom, 
we are led to look for an eigenstate of 7i in the form 

tl> {1) (si) = (^j e liai+bl)si (Aie iqSl + Bie~ iqSl ) 
+ (° j e l(ai - bl)si (Cie iqSl +Die~ igSl ), 
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V> (2) (s 2 ) = ( C ^f)e^ a2+b2)s2 (A 2 e iqS2 + B 2 e~ iqS2 ) 

&*- b2 >2(C 2 e lqS2 + D 2 e- lqS2 ). (8.6) 

Along the lines of Section [2J the continuity conditions generalizing (12.51) allow one to 
express the eight amplitudes Ai, . . . , D 2 in terms of the two components of ^(P). The 
current conservation conditions generalizing (12.61) then yield the characteristic equation 
D$o{q) = 0, with 

Dso(q) = (sing(Li + L 2 ) - cos($ 2 + ^ 2 ) singLi - cos($i + \I>i) singL 2 
x \sm.q(L\ + L 2 ) — cos($ 2 ~ ^2) singi^ — cos($i — ^i) singL 2 

+ 4 singLi singL 2 sin$i sm$ 2 sin^i sin\l/ 2 sin 2 (8.7) 

The first two lines of this expression are identical to the scalar characteristic 
equation ( 12.1 Oft , up to the replacement of the magnetic fluxes by the sums and differences 
of their Abelian and non-Abelian parts: $1 — > $1 ± ^1, $ 2 — > $ 2 ± \l/ 2 . Each factor 
therefore describes a scalar problem in an effective Abelian flux. The third line provides 
the coupling between both spin components. In order for this coupling to be non-zero, 
the four fluxes need to be simultaneously non-trivial, i.e., not equal to or ir (mod. 2ir), 
and the angle j3 not equal to (mod. 2%). 

In the case of two equal rings (L 1 = L 2 = L, $1 = $ 2 = $, ^ = ^ 2 = the 
spectrum consists of an alternation of groups of two degenerate idle states, such that 
singL = 0, i.e., qL = pir, and of groups of two non-degenerate active states, such that 

cos qL = cos9± = cos $ cos \I/ ± sin $ sin \1/ cos ^, (8-8) 

Choosing for definiteness < 3 < n, we have < 9 + < 9- < ir. For N = 2 electrons, 
only the first two active states are occupied. The magnetization reads 

M = -LL°4± + e JlA. (8 . 9) 



<9$ / 

Let us consider the dependence of M on the Abelian flux $ in the range < $ < 7r 
and for < ^ < 7r. For (3 = 0, we obtain 

M =(7* / " . ° (8.10) 

( (7T - $)/% for 7T - ^ < $ < 7T. 

This discontinuous jump in the magnetization is rounded for non-zero values of /3, as 
shown in Figure [HJ This figure also illustrates a remarkable feature of such a simple 
system of two electrons on two equal rings, namely that the magnetization changes sign 
as a function of parameters. For a small magnetic flux $, we have 

M = — — (cos 2 ^ + sin 2 f ¥cot¥) $ + •■■ (8.11) 
7T \ 2 2 J 

The expression in the parentheses vanishes for 
9 8 tan \1/ 

tan 2 - = — , 8.12 

2 ^ ( J 
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$/rr $/tt 

Figure 8. Plots of the magnetization M of two electrons on two equal rings, against 
<&/7r, for * = 7r/2 (left) and * = 7r/5 (right). In both cases, /3 = kir/10, with 
k = 0, . . . , 10, bottom to top in the left part of the curves. 



provided \& > tt/2. More generally, the magnetization vanishes along a ^/-dependent 
curve in the ($,/3) plane, shown in Figure [S At fixed the magnetization is positive 
above and to the right of that curve, whereas it is negative below and to the left of that 
curve. 




0.2 0.4 0.6 0.8 1.0 

$/tt 



Figure 9. Plot of the locus in the (<£>, /3) plane where the magnetization M of two 
electrons on two equal rings vanishes, at fixed \&. From top right to bottom left, 
*/tt = 0.1, 0.2, 0.3, 0.4, 0.45, 0.49, 0.5 (vertical line), 0.51, 0.55, 0.6, 0.7, 0.8 and 0.9. 



9. Discussion 

The many differences between the single and the double-ring geometry, which have 
been underlined throughout this work, testify the importance of topology in mesoscopic 
physics. From a local viewpoint, both systems are governed by the same one- dimensional 
Schrodinger equation. They only differ in their global topological structure, measured 
by their genus (number of holes) . 
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Another example underlining the key role of topology in the context of persistent 
currents has been studied by a collaboration involving one of us [15] . The importance of 
the sample's topology in mesoscopic physics has already been stressed by Schmeltzer [TO] , 
who studied the persistent current of spinless fermions in a double-ring system, using the 
formalism of Dirac constraints. In fact, Schmeltzer's work was the main motivation for 
the present research. Apart from underlining the importance of topology, the work [TO] 
however focussed on different aspects. In particular, the role of spin and the AC effect 
were not tackled. 

We have not tackled the rich physics that should emerge when interaction effects 
are taken into account. First of all, there are interesting charging effects, suggested and 
discussed for the single- ring geometry in [161 E] • Second, interactions in one dimension 
turn the system to a non-Fermi liquid. While this problem has been studied for a single 
ring [18], persistent currents in interacting quantum rings have recently been studied 
in [19]. It can however be anticipated that the detailed analysis presented here for the 
non-interacting system will turn its place to a complicated and somewhat intractable 
formalism, so that e.g. the fine effects of length and flux commensurability will be absent. 

The two-ring geometry provides a playground for testing fundamental aspects of 
quantum mechanics, such as the occurrence of interlaced AB and AC effects. At the 
same time, a two-ring system is certainly within the reach of fabrication (see Figure 1 
in [20]), so that the present study is also rooted into the real world. At the same 
time, the difficulties of controlling the external electric field in semiconductor he- 
terostructures hosting a two-dimensional electron gas, as well as their enhancement 
by internal fields of ion cores and discontinuities (eventually manifest as spin-orbit 
couplings) have recently been summarized in a comprehensive review article [2TJ. Thus, 
while the bare Hamiltonian (18.41) of course applies for any field whatsoever, a realization 
of the two-ring geometry as in Figure[7]is still to be considered as a Gedankenexperiment. 
Nevertheless, the main conclusion of our analysis is fundamental and should pass any 
experimental test, namely, an AC effect which is a pure gauge (periodic in the SU(2) 
phase) and does not conserve s z is realizable only if the AB effect is present. 

Finally, although the present work has focussed onto zero-temperature equilibrium 
properties, it can be expected that transport properties will reveal a similar richness of 
behavior. 
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Appendix A. One single ring: a reminder 

In this Appendix we give a brief reminder of the well-known situation of a single clean 
ring. The ring with length L and area A is threaded by a flux $ = BA. It may assume 
an arbitrary planar shape. With the same conventions as in the body of this work, 
parametrizing a point of the wire by its curvilinear abscissa s in the range < s < L, 
the one-body Hamiltonian reads 

H = (p-a) 2 , (A.l) 

with p = — ld/ds and a = 

It is worth mentioning the analogy with Bloch theory, already noticed in [2]. 
Starting from the Schrodinger equation Tiip = Etp with ip(s + L) = ip(s), a gauge 
transformation ijj(s) = e ms r](s) leads to the following equation and boundary condition 
for rj(s): 

p 2 r](s) = Erj(s), 7]{s + L)=e 1 *r]{s). (A.2) 

This is exactly the equation for an electron in a one-dimensional lattice potential of 
period L and Bloch wavenumber K such that KL = $, i.e., K = a = $/L. One 
immediate consequence is that the energy is periodic in $ with period 2ir. This results 
also holds when the ring is not clean, since the full (disordered) potential can still be 
viewed periodic potential of period L. 

The eigenstates of H are given by the wavefunctions 

t/}(s) ~e l( " +a)s . (A.3) 

The periodicity of tp(s) in s with period L yields the quantization condition (q + a)L = 
2trk, hence 

Qk = ^ , (A.4) 

where k = 0, ±1, . . . The wavefunctions ipk{s) ~ exp(27ri/cs) therefore do not depend on 
the magnetic flux $. 

The contributions and of the eigenstate number k to the persistent current 

and the magnetization read 

dE k 2q k dE k t . . 

M k = -— ± = AI k . (A.5) 



9$ L ' dB 
With the notation (I3.3p . the above result for I k simply reads Q k = 1. 

For a zero-temperature system with electrons, and for < $ < ir, we have the 

following results. 

• For N = 2p — 1 odd, the occupied states are k = —p + 1, . . . ,p — 1. We obtain 

N 2A 
E = E + — $ 2 , M = -— N$. (A.6) 

• For N = 2p even, the occupied states are k = —p + 1, . . . ,p. We obtain 

AT OA 

E = E + — (7r-^) 2 , M=— A(tt-$). (A.7) 

In both cases the minimum energy reads E = N(N 2 — l)ir 2 /(3L). For a circular 
ring with radius R, we have L = 2ttR and A = irR 2 , so that 2A/L 2 = l/(27r). 
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Appendix B. Three coupled rings 

In this Appendix we show how the present investigation can be extended to more 
complex geometries. We consider for defmiteness the case of a sample made of three 
unequal rings touching at two contact points P and Q, as shown in Figure EU The line 
PQ is assumed to be an axis of symmetry of the sample. The left, middle and right rings 
have respective lengths L%, L 2 and L 3 and areas A%, A 2 and A3. They are therefore 
threaded by magnetic fluxes $1 = BA\, $2 = BA 2 and $3 = BA 3 . 




Figure Bl. The sample considered in Appendix B is made of three unequal rings 
touching at two contact points P and Q. 



The one-electron Hamiltonian of the system reads 

H = {pi - ai) 2 + ip 2 - a 2 f + (p 3 - a 3 ) 2 , (B.l) 

withpi = -id/dsi,p 2 = -id/ds 2 ,f>3 = -id/ds 3 , a x = $i/Li, a 2 = $ 2 /L 2 , a 3 = $ 3 /L 3 . 

A state is now described by three wavefunctions, one living on each ring: 
{ip^(si), ip( 2 >(s 2 ), ip^(s 3 )}. These wavefunctions obey two continuity conditions of 
the form (12.51) (one at P and one at Q) and two current conservation conditions of the 
form (12. 6p . Along the lines of Section [2], we obtain the characteristic function 

D 3 (q) = cosg(Li + L 2 + L 3 ) 

— cos $ 3 cos q(Li + L 2 ) — cos $1 cos q(L 2 + L 3 ) 

+ cos $ 3 cos qLi + cos $1 cos $ 3 cos qL 2 + cos $1 cos qL 3 

+ cos $2 sin qL\ sin qL 3 — cos qL\ cos qL 3 — cos $1 cos $ 3 . (B.2) 

Many of the outcomes of this paper can be extended to the present situation, 
although expressions become very cumbersome. The energy eigenvalues E n = q\ can be 
parametrized as 

Qn = T T T (n = 1,2,...), (B.3) 
J^l ~r Li 2 -+- ±j 3 

where the modulation g n obeys the bounds — ir < g n < 2ir. In the absence of magnetic 
fluxes, the factorized form (14.1j) of the characteristic function generalizes to 

n / n a . g(^i + ^ 2 + L 3 ) . qLi . qL 2 . qL 3 

L> 3 (q) = 8 sm sm sm — — sm . (B.4J 

The spectrum consists of four sectors. The corresponding states can be respectively 
referred to as trilateral, left, central and right. 
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